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We explore the motion of a classical particle in a symmet- 
ric potential with non-Gaussian skewed white noise. We show 
analytically and numerically that the presence of nonzero odd 
moments leads to a macroscopic current. For a noise with a 
vanishing third moment we find that the current changes sign 
as the potential or the noise strength is increased. Possible 
physical situations are discussed including motor-protein mo- 
tion and driven systems like fluidized beds. 
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There can be no current in a system in thermal equi- 
librium. Drift is usually associated with the presence 
of an external field, for example an electric field as in 
the case of charge drift. However, it has been shown 
in recent studies 0] that a spatially asymmetric poten- 
tial (usually taken to be an asymetric sawtooth) which 
exerts no macrosopic force, combined with temporally 
periodic forcing or a noise of non-thermal origin leads 
to non-equilibrium steady states that carry current. It 
is believed that these models may have relevance in bio- 
logical transduction and the motion of motor proteins in 
eukaryotic cells. 

In this paper we show that (1) it is not necessary to 
have a spatial potential that violates parity to get uni- 
directional drift [^|-^|; a combination of symmetric po- 
tential and non-thermal skew noise is sufficient; (2) the 
current vanishes in the limits of very strong or very weak 
potentials or noise amplitudes (Fig. [j]); (3) if the skew 
noise is dominated by a third moment then the current 
always has the same sign as the third moment (Fig. [I], 
inset); (4) higher moments can lead to current reversal 
(Fig. ^). Finally we construct a Maxwell demon based 
on our model and also suggest an experiment involving 
a colloidal particle in a fluidized bed. 

Perhaps the simplest physical example of a skew noise 
occurs in the case of a Brownian particle falling under 
gravity in an ideal gas. The noise is the force in the rest 
frame of the particle due to random kicks coming from 
collisions with gas molecules. Since the net force in the 
rest frame is zero, the kicks in the downward direction, 
though fewer in number, have to be stronger, leading to 
a skew noise. Besides this specific example, the noise is 
expected to be skewed in any driven system where the 
driving violates parity, e.g. colloidal suspensions in a 
fluidized bed, driven diffusive lattice gases etc. 

To capture the essentials of the problem, we consider 
overdamped motion of a particle in a periodic symmetric 



sawtooth potential V(x) with a skewed noise. Suppose 
the noise kicks come at regular intervals of time A. The 
displacement in one such interval, Sx = a + 77, has a de- 
terministic part a — -T^A where T is the inverse of 
the friction coefficient, and a stochastic part 77, defined 
by: 



r\ = —s/3 with probability pi 

7/ = +2s/3 with probability p r 

7/ = with probability 1 — p m 

The noise is uncorrelated in time by construction. p — 
(1 — p m ) is the probability that the particle does not ex- 
perience a noise kick during the interval. 
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FIG. 1. Variation of current with strength of confining po- 
tential: Monte Carlo Simulation in symmetric sawtooth po- 
tential. The different curves correspond to different values of 
the noise strength. The probablity of moving to the left and 
right respectively in each step was pi = 0.2 and p r = 0.1. The 
current plotted is the distance moved in one timestep, where 
the period of the potential is = 1. The inset shows the current 
in the approximate analytic solution of the master equation 
for the same set of parameters. 

If the stochastic part of the displacement is much 
smaller than the the length ly over which the determin- 
istic force changes appreciably (ly = 1 in our case), then 
it is legitimate to consider a larger timestep TV A, within 
which N noise kicks are present, adding up to a Gaus- 
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sian thermal distribution which can cause no drift. If, 
however, the displacement due to noise is comparable 
to ly, the local environment changes considerably after 
each noise-kick, so that it is no longer permissible to add 
several realizations and the non-Gaussian nature of the 
noise manifests itself. In the extreme limit where the po- 
tential is weak and the noise strength s lies between i 
and 1 (where the potential has a period 1), it is easy to 
see qualitatively that there will be a drift in the direc- 
tion of the stronger but less frequent kicks. It turns out 
that even for the less obvious case of weaker noise am- 
plitudes, there is a current towards the right, which goes 
continuously to zero as the noise amplitude is decreased. 

Direct Monte Carlo simulations have been carried out 
on this model (Fig. [j]). The simulation was run for a 
long enough time (10 8 steps) so that the drift (which 
goes linearly in time) was significant compared to the free 
diffusion (i.e., the sum of the random numbers), which 
goes as y/t, and the current had acquired a steady value. 
The drift current is expected to go to zero in two limiting 
cases for different reasons. When the potential is weak or 
the noise strength is large, the system behaves as though 
there were no potential. On the other hand, when the 
potential is strong or the noise is weak, the particle is 
confined. The current is always positive, i.e., along the 
direction of the stronger kicks. This, as we shall see be- 
low, is a result of the third moment of the noise playing 
the dominant role. 

We now look at the evolution of the probability density 
in the above problem. First consider a general Markoff 
process, one in which the random noise can have mo- 
ments higher than the second. Let P(x,n)dx be the 
probability of finding the particle between x and x + dx 
at timestep n. The time evolution of P(x, n) is given 
by P(x, n + 1) = J duP(x — u, n)T u (x — it), where T u [x) 
is the probability of transition from x to x + u in one 
timestep. The functions inside the integral on the right 
hand side can be Taylor expanded about x to give an 
evolution equation for P involving the moments of a; [|| . 
The rate of change of probability is the divergence of a 
current, given by P(x,n + 1) — P(x,n) = — -§zJ, where 
the current is 



J={u)P(x,n)-^[(u 2 )P(x,n)\ 



1 d 2 



[(u 3 )P(x,n)]. 



(1) 



The moments (u m ) (x) = J duu m T u (x) that enter into 
this equation are the combined effect of the potential 
and the noise. For the model described above, where 
the displacement in timestep A is a(x) — — §j-A; due to 
the potential and — Si and s r with probabilities pi and p r 
respectively due to the noise, the moments are (u) = a; 
(u 2 ) = a 2 + fi2'- (u 3 ) = a 3 + 3a/i2 + P-3', and so on, where 
H n is the rfi 1 moment of the noise: \i n = pi(—si) n +p r s™. 



For an arbitrary potential, a is a function of x whereas 
the fi n s are constants. 

For ordinary Brownian motion, the series in equation 
([j]) may be truncated at the second term. To see that 
there is no current in this case, simply divide equation [j] 
(truncated at the second order) by P and integrate both 
sides from to 1, with the limits identified with each 
other due to periodic boundary conditions. The first mo- 
ment is caused by the static potential alone, since noise 
has a zero mean. It is antisymmetric about ^ for a sym- 
metric potential, so its integral is zero. The second term 
is a total derivative, therefore it's integral also vanishes 
on a circle. Since the integral of -p on the left hand side 
is positive definite, J must be zero. 

Let us now restrict ourselves to processes in which the 
moment is given by (u n ) = ji n . This amounts to 
assuming that a/s is small, so that the entire contribu- 
tion to the second and higher moments comes from the 
noise alone, an approximation which holds good for suf- 
ficiently weak potentials. The current can can again be 
calculated as above, but the additional terms on the RHS 
are no longer total derivatives. We also assume that the 
higher moments p n of the noise are negligible compared 
to the third moment /Z3 . Under these approximations the 
current turns out to be 
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which clearly has the same sign as the third moment /X3 
in agreement with intuitive expectations as well as our 
numerical results. 

When the third moment is much smaller than the sec- 
ond, i.e., for slight deviations from a thermal distrbu- 
tion, we can calculate this current perturbatively. The 
probability distribution is the equilibrium distribution 
plus a correction due to the third moment, P(x) ~ 
e -pv(x) _|_ (fj, 3 ). Since the current is already first order 
in /i 3 , we can use the equilibrium probability distribution 
itself in equation (0) to get 
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As expected, this current goes to zero when the potential 
becomes flat. Curiously, if we add a slight disorder to the 
potential, thereby making it more corrugated, i.e., give 
it a larger value of (^r) we get an enhancement in the 
current. 

To express the current (|3|) in terms of the parame- 
ters a and s, we compare the small-kick limit of the 
discrete model with the equation of motion of a parti- 
cle in a potential V and temperature T. This gives us 
the connection: a\x\ = TAV(x) for the potential and 

2Tk B TA = (f^-) for the noise. Equation (g) then 
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This cur- 



gives J = ^ 3 A 3 /(e A - 1) where A = 
rent is shown in Fig. [l] (inset). In spite of the drastic 
approximations there is excellent qualitative agreement 
with the numerical result. As expected, the drift van- 
ishes in the limiting cases of (i) weak potentials or strong 
noise and (ii) strong potentials or weak noise, in a manner 
which we now calculate. 

For A — ► and we have J ~ /^(pr) showing 
that J ~ a 2 in the limit of weak potential strength. 
Since, ^ ~ s 3 the current falls off as s _1 as the noise 
strength s is made large. On the other hand, for large 
A, J ~ jjexp _ constant _ Tnig ig tng manner in w hi c h 

the current approaches zero as the noise kicks s become 
smaller and the validity of the central limit theorem im- 
proves. 
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FIG. 2. Reversal of the sign of current as potential strength 
is increased: Monte Carlo simulation on symmetric sawtooth 
potential with a noise with three kicks such that first and 
third moments vanish. The noise strength s for this run was 
0.1 and the noise was composed of three kicks of strength 
s, —2s and 3s with probabilities \p m , \Pm and ^p m , with 
p m = 0.3. The run was for 10 8 steps. 

Interesting features develop when the noise is domi- 
nated by the fifth moment. If we neglect all other mo- 
ments, The current can again be calculated perturba- 
tively: 
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/? 2 (0) ; 
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The second term would lead to a current opposite in sign 
with respect to the fifth moment. For a sawtooth po- 



tential, the integral of (j^r^J is ill defined. If, however 
we replace the sawtooth by a potential with a rounded 
wedge so that the derivative ^ changes from —a to +a 
over a length e, then the current is 



J = 
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Evidently, the negative part of the current dominates for 
small A, i.e., for weak potentials a or large noise strength 
s. The sign changes as the strength a is increased. We 
have found such a reversal in direct simulations as well 
( Fig. |2 _). Current reversal has been found in the lit- 
erature |9|] for a particle in an ratchetlike potential with 
coloured noise. In the present model, the presence of 
higher moments in a white noise causes current reversal 
as a parameter, e.g., the potential strength is changed 
continuously. 

Having explored the consequences of the model in some 
detail, we now go back to physical realizations. Before 
going on to a concrete experimental situation where skew 
noise may arise, let us perform a thought experiment. 
Consider two 'gases' having different particle masses kept 
in separate containers, say Hydrogen (lighter atoms) and 
Nitrogen (heavier). Let both the containers be at the 
same temperature and pressure, and hence the same den- 
sity (the densities are sufficiently low so that the gases 
can be treated as ideal). Now connect the reservoirs by 
a long tube which has a partition in the middle sepa- 
rating the two. The pressure is the same on both sides 
so the partition remains stationary. Note, however, that 
the random force that the partition experiences is a noise 
precisely of the kind discussed above JT3] . Now, if there 
are minute corrugations on the base of the tube then they 
supply a periodic potential due to the gravitational field, 
and the partition should start moving into the lighter 
gas. It is possible even to tilt the tube slightly and drive 
the partition up against the gravitational field, thereby 
extracting work from a single system at thermal equilib- 
rium! 

The demon can be rationalized along lines similar to 
other more familiar ones [jjLlf] . We have assumed over- 
damped dynamics of the partition and neglected the noise 
associated with it, which is always there in any source of 
dissipation. If the corrugated floor is colder than the 
reservoirs, then this noise amplitude ( which, for a given 
damping, is proportional to the square root of tempera- 
ture) is negligible, and the partition drifts. In this sit- 
uation we are driving an engine between the gases at a 
higher temperature and the floor at a lower temperature. 
As the partition moves, it keeps transferring heat to the 
floor, until it reaches the temperature of the reservoirs. 
Now the thermal fluctuations of the partition completely 
wash out the effect of the asymmetric noise (coming from 
collisions) and the drift stops. 

A more concrete experimental situation may be real- 
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ized using charged colloidal particles in a fluidized bed. 
A steady upward flow of fluid supplies the Stokes friction 
that is adjusted to balance exactly the weight of the parti- 
cle. Brownian motion is negligible for most experimental 
systems, but there are indications of the existence of a 
'noise' on the particle due to hydrodynamic interactions 
p2[ . The source of this noise is completely athermal. 
From symmetry considerations, this noise is expected to 
be skewed, in a way similar to the noise on a falling parti- 
cle in an ideal gas. If we now switch on a sinusoidal laser 
field forming standing waves with an intensity variation 
in the vertical direction, the particles will start drifting 
downwards in the lab frame. The drift speed will depend 
on the flow rate in the fluidized bed (which controls the 
noise amplitude s) and the amplitude of the laser field 
(~a). 

Finally, a word about biology. The 'standard model' 
with a ratchet-like potential and some kind of correlated 
noise has been advertised as a candidate for modelling 
motor protein motion in eukaryotic cells. However, there 
appears to be no direct evidence of such ratchet poten- 
tials in actual experiments. The model presented here 
can be an equally well suited candidate for motion of 
proteins and has all the essential features of the stan- 
dard model, including a possible reversal of current as 
a parameter is continuously changed. The question of 
current reversal is interesting because motor proteins of 
the same class are known to move along different direc- 
tions; for example a drosophilla kinesin, Ned, and axonal 
kinesin walk in opposite senses along a microtubule [jl3| . 
Our model gives, for the case where the fifth moment 
dominates, a reversal of current simply by changing the 
noise amplitude (which should correspond to ATP activ- 
ity) leading to the possibility of two-way traffic along the 
highways. 

In conclusion, we have shown that it is possible to pro- 
duce a macroscopic current in a parity non- violating pe- 
riodic potential if the noise has nonzero odd moments. 
The current vanishes when either the potential strength 
or the noise amplitude become very small or very large 
and exhibits sign reversal if the third moment is negli- 
gible. The vanishing of the current in the strong noise 
limit distinguishes our model from that of Luczka et. al. 
Pj who consider a symmetric potential and a shot noise 
and find a saturation of the current in the strong-noise 
limit. This is because of the extreme nature of the shot 
noise such as the divergence of moments higher than the 
second. 
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